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ABSTRACT. A codimension-three bifurcation, characterized by a pair of purely 
imaginary eigenvalues and a nonsemisimple double zero eigenvalue, arises in 
the study of a pair of weakly coupled nonlinear oscillators with Z2 © ^2 sym- 
metry. The methodology is based on Arnold's ideas of versal deformations of 
matrices for the linear analysis, and Poincare normal forms for the nonlinear 
analysis of the system. The stratified subvariety of primary bifurcations of 
codimensions one and two is identified in the parameter space. The analysis 
reveals different types of solutions in the state space, including equilibria, limit 
cycles, invariant tori and the possibility of homoclinic chaos. A mechanism is 
identified for energy transfer without strong resonance between two oscillation 
modes with widely separated frequencies. 



This work is motivated by the study of weakly-coupled nonlinear oscillators 
which may be modeled by second-order differential equations of the general form 



Here j\ and f 2 are holomorphic functions of the state variables x, x, y, y and param- 
eters n G The functions are higher order in that they contain only quadratic 
and higher order terms in these state variables and parameters, which are assumed 
to be small, and fj(0, 0, 0, 0, /i) = 0, j = 1,2. The parameters /j, represent physi- 
cally relevant quantities such as damping, detuning and the coupling between the 
oscillators; the role of these parameters will be clarified further below. In (1.1), the 
natural frequencies of the two oscillators, when linearized at (x, x, y, y) = (0, 0, 0, 0), 
are given by W12 > 0. Systems of the form (1.1) occur frequently in the study of 
mechanical systems, electrical circuits, the biological sciences and elsewhere. 

A major concern in the study of coupled oscillators is the understanding and 
control of resonance phenomena. V.I. Arnold has made influential contributions to 
the general theory of resonances, for example [5, 8, 9]. The ratio of the natural 
frequencies plays a leading role in determining resonances in systems (1.1) of two 
oscillators. The cases = -, where p and q are small positive integers, p + 
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1. Introduction 




x + u\x + fi (x, x, y, y,fi) =0 
y + u)'iy + f 2 {x,x,y,y,ii) = 0. 
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q < 5, are called strong resonances. The remaining rational cases are called weak 
resonances, and cases with the ratio — irrational are called nonresonant. As the 
names suggest, coupled oscillators interact strongly at or near strong resonances, 
and interact weakly, if at all, in the other cases. This is explained in part by the 
fact that in the Poincare normal form corresponding to these coupled oscillator 
equations, terms of degree p + q in the state variables remain. These are called 
resonant terms. Locally, these terms have a strong influence when p + q is small 
and their influence decreases as p + q increases. 

This paper addresses the case of coupled oscillators (1.1) with natural frequen- 
cies which are widely separated; we assume u>i very small and u>2 very large, relative 
to each other. On the basis of the above classification, this would not be considered 
a case of strong resonance. To fix magnitudes, rescale time in (1.1) so that u>2 = 1- 
Then our assumption is u>i — > 0. Equivalently, we consider natural frequencies at 
or near the ratio — = t- Therefore, we refer to this as the case of : 1 resonance. 

In this paper, system (1.1) is viewed as a perturbation of a singular system with 
a bifurcation of codimension three, arising through the coalescence of two simpler 
bifurcations. The low frequency mode, in the limit u>i — > 0, is a Bogdanov-Takens 
bifurcation [11, 41]. This bifurcation has codimension two. The high frequency 
mode is close to a nondegenerate Andronov-Hopf bifurcation, which has codimen- 
sion one. Historically, this bifurcation has its origins in the work of Poincare [37, 
sections 51-52,316] more than a century ago, and was studied extensively by An- 
dronov and coworkers [1, 2, 3, 4] starting in the 1920's. Arnold [9, section 33. A, B] 
presents a pellucid account of the Poincare-Andronov theory. An excellent account 
of the approach of Hopf and later developments is given in [30]. Degenerate cases 
of Hopf 's theorem are explored in [17]. Although many different names have been 
used for this theory over the years, in this paper we choose the name Andronov- 
Hopf bifurcation, which is suggestive of its rich history. 

Using Arnold's theory of versal deformations of matrices [6], we verify that 
the singular system combining these two classical bifurcations has linear part with 
codimension three. Then we investigate via Poincare normal forms the implications 
of this codimension-three singularity for nonlinear systems. The analysis reveals, 
in addition to features associated with the original two classical bifurcations, new 
phenomena arising from the nonlinear mode-interactions. 

Cases of bifurcations involving a simple zero eigenvalue with a pair of purely 
imaginary eigenvalues (and no other eigenvalues with zero real part) have been 
studied by Arnold [9], Langford [24, 25], Guckenheimer and Holmes [18], Khazin 
and Shnol [23], Chow et al. [14], Scheurle and Marsden [39] and others cited in these 
references. Less is known about the case of a nonsemisimple double zero eigenvalue 
together with an imaginary pair. Some early analytical results and applications 
have been presented by Moson [31, 32], Yu and Huseyin [45], Georgiou et al. [15], 
Nagata and Namachchivaya [33], and Nayfeh et al. [34, 35, 36]. 

In a series of papers [34, 35, 36], Nayfeh and coworkers have investigated a 
system of two weakly coupled oscillators with widely separated frequencies similar to 
(1.1). They point out an interesting paradox. In both mechanical experiments and 
a mathematical model, when the system is excited near a high natural frequency, 
large low-frequency responses accompanied by slow modulations of the amplitude 
and phase of the high-frequency mode are observed. Thus, they found that there is 
a strong interaction between these two modes, and a transfer of energy from high 
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to low frequency modes, even though by the above classification it is not a case of 
strong resonance. Their work provided additional motivation for the present paper. 

To fix ideas, throughout this paper it is assumed that each of the oscillators 
in (1.1) has an odd reflectional symmetry about the rest position. We extend this 
symmetry to the coupled system by assuming that the system commutes with the 
TL'i © TL'i symmetry group generated by two reflection operators 

f x {-x,-x,y,y,n) = -f±(x,x,y,y,n) 
f 2 {-x,-x,y,y,n) = f 2 {x,x,y,y,n) 



(1.2) 



fi(x,x,-y,-y,/i) = h{x,x,y,y,n) 
f 2 {x,x,-y,-y,n) = -f 2 (x,x,y,y,n) 



This symmetry is very common in applications; for example, the classical pendulum 
equation, Duffing 's equation and Van der Pol's equation all have odd symmetry, as 
do the examples in [33, 34, 35, 36]. The symmetry (1.2) forbids linear coupling 
terms (i.e. linear x, x terms in the {/-equation or vice-versa) but permits nonlinear 
coupling (such as x 2 y in the {/-equation). Systems with less symmetry that (1.2) 
may also be studied by the methods of this paper. 

Before beginning to analyze (1.1), we rewrite it with the linear part in a more 
explicit form. Equation (1.3) displays the most general linear terms in the state 
variables, which are consistent with the symmetries (1.2) and the time-rescaling. 
(Recall that t has been rescaled so that u> 2 = 1). The functions fj, which replace 
the former functions fj , consist of all the remaining nonlinear terms in the state 
variables (x, x, y, y), with coefficients which are functions of the parameters \i. Since 
the functions fj are assumed to commute with the two reflection operators (1.2), 
they are cubic or higher order in the state variables. 

^ 3 -j x + Six + Clx + fi(x,x,y,y,fi) = 

jj + hy + y + f 2 {x,x,y, y, fi) = 

In the linear terms of (1.3), S\, S 2 are linear damping parameters and £1 = (oji/oj 2 ) 2 
is a small tuning parameter which from here on is allowed to be negative or zero; 
therefore, £1 is better thought of as a restoring force than a frequency ratio. In the 
engineering literature, a change in sign of SI from positive to negative is associated 
with a transition from "flutter" to "divergence" . These physical parameters are all 
assumed to be small. Recognizing the physical importance of these three parameters 
and choosing to ignore the possible effects of other parameters, we formally define 
the parameter vector /j, to have precisely these three components 

(i.4) fi = (n,s 1; s 2 ) . 

The appropriateness of this choice will become clearer after the presentation of 
Arnold's versal deformation, in the next section. 
With the definition 

f yi \ I x \ 
2/2 _ x 

ys ~ y 

V f4 / V y ) 

system (1.3) can be written as the following system of first order equations, with 
linear part which is a perturbation of a matrix Ao in real Jordan form. 

(1-5) y= [Ao + CMJy + f , 
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Here C(ji) is the linear perturbation matrix, depending only on the three small 
physical parameters \i = (£1, Si, S 2 ). 

We may simplify notation by introducing a mixed set of real/complex state 
variables with (t/i, y 2 ) £ IR 2 defined as above, and z £ C defined by 

(1-6) z = y 3 - iy 4 , z = y 3 + iy 4 . 

(The reason for the nonstandard choice of signs of the imaginary parts of z, z is the 
following. In the phase plane of a second order harmonic oscillator equation, the 
coordinate axes are normally chosen to be (x,y) = (x,x), which forces oscillations 
about the origin to have a clockwise orientation. In the complex plane, the conven- 
tion is that a counterclockwise rotation is positive, as for the differential equation 
z = iz. The transformation (1.6) effects this reversal of orientation, between the 
real and complex phase planes.) Then system (1.5) is transformed to the system 

2/1 = 2/2 

(1-7) y-i = -fii/i - Siy 2 - fi 

z = (i - \&2)z + \S 2 z + if 2 

together with a fourth equation for z obtained by conjugation of the equation for 
z given here. 

The matrix Ao in (1.5) is recognized as being in real Jordan form, with the 
four eigenvalues 2 (nonsemisimple) and ±i (imaginary pair). The correspond- 
ing differential equations (1.5) for small amplitudes represent a perturbation of a 
system of two independent oscillators, one with coordinates (x, x) and natural fre- 
quency u>i 0, and the other with coordinates (y, y) and natural frequency u> 2 = 1. 
However, the results of this paper apply more generally than this. Consider any 
parametrized two degree-of- freedom system, which may have coefficient matrix far 
from Ao; but suppose that, for some value of the parameters, the linearization at 
an equilibrium has eigenvalues {0 2 , ±i} (after rescaling time). Then the linear part 
can be transformed to Ao by an invertible linear change of coordinates (similar- 
ity transformation), which may be far from the identity. In this case, the state 
variables in (1.5) are the normal mode coordinates of the original system, and the 
symmetries (1.2) are symmetries of the normal modes. The state variables in this 
paper should be understood in that light. 

Even more generally, the results of this paper apply to systems of ordinary 
differential equations of dimension higher than four, and to classes of infinite di- 
mensional partial differential equations or functional differential equations, satisfy- 
ing appropriate spectral conditions. If the linearization at an equilibrium solution 
satisfies standard compactness and hyperbolicity assumptions on the linear opera- 
tor and its spectrum, then the Center Manifold Theorem applies. In this case, if 
the spectrum contains the eigenvalues {0 2 ,±i}, and the complementary spectrum 
is bounded away from the imaginary axis in the left half-plane, then the study 
of long-time dynamics of the full system reduces to the four dimensional system 
considered here. Such higher dimensional generalizations are not pursued further. 
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The paper is organized as follows. In the next section, a versal deformation 
of the linear part of the singular vector field is obtained, using ideas of Arnold, 
and the codimension is computed to be three. Section 3 presents the Poincare 
normal form for this system. The original four-dimensional system is reduced to 
a three-dimensional one, exploiting the S 1 normal form symmetry arising from the 
Andronov-Hopf bifurcation. In Section 4, the generic bifurcations of codimension 
one and two are explored. A rich variety of solutions, including stationary, periodic 
and quasiperiodic (tori) are found. The main results of this paper include the strat- 
ified subvariety of primary bifurcations, and formulas for primary and secondary 
bifurcations of equilibria, periodic solutions and tori. In Section 5, the relevance 
of these results to the energy-transfer phenomenon of Nayfeh [34, 35, 36] is con- 
sidered by means of an example, and directions of further research are indicated in 
Section 6. 

2. Linear Analysis 
Definition. Let C™ denote the set of n x n complex matrices. A family A of 

2 

matrices on the base M is a holomorphic mapping A : M — > C™ , where M is a 
neighbourhood of the origin in the parameter space C fe . For brevity, we call this 
the family A. The germ of a family A at the point 6 M is called a deformation 
of the matrix A(0). The elements /j, £ M are called parameters. 

An example of a family of real matrices is given by Ao + C(jJ-) in equation (1.5) 
and a complex family is given by the linear part of (1.7). The algebra is simpler 
if we work with complex matrices. The real Jordan matrix Ao of (1.5) has the 
complex Jordan form 

10 \ 


i 
-i ) 

Two main questions are addressed in this section. The first question is how 
to find a canonical matrix, which is as simple as possible, and to which the linear 
part of all coupled-oscillator systems as described in Section 1 can be reduced by 
a linear change of coordinates. An obvious candidate is the Jordan form; however, 
this is not a good choice because the reduction of a matrix to its Jordan form 
is not a stable operation when there are multiple eigenvalues, as is the case here. 
Although every member of a family containing A c may be reduced to a Jordan form, 
the Jordan form does not in general depend continuously on the parameters, since 
there are matrices arbitrarily close to A c for which the Jordan form is diagonal. The 
mapping which transforms matrices to Jordan form is discontinuous. Therefore, we 
seek a different normal form, to which members of any family containing A c may 
be transformed smoothly. 

The second question relates to the choice of parameters /j, £ M. What is 
the minimum set of parameters needed, and how can we know if any essential 
parameters are missing? In particular, are there systems with linear parts close to 
Ao which are not attainable by variations of the physical parameters in the given 
perturbation matrix C(ji) in (1.5)? If we have in fact overlooked some essential 
parameters, then any study of the system in Section 1 is likely to miss important 
phenomena. We would like to have assurance that the parameter set in our family 
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of matrices is sufficiently rich that all perturbations of Ao are included, and that 
there are no redundant parameters. On the other hand, it may happen that due 
to the mechanical structure of equations (1.3), the parameters and forces can enter 
only in certain ways and generic perturbations are unrealistic for the system. In 
that case, a comparison with the generic case is useful, to uncover any previously 
hidden constraints in our model. 

Both of these classes of questions were answered by Arnold [6, 7, 9], who 
introduced a generalization of the Jordan form of a matrix A, with the minimum 
number of essential parameters, and to which all members of any family of A can 
be transformed smoothly. This generalization of the Jordan normal form may 
be called the Arnold normal form of a matrix. The remainder of this section 
summarizes Arnold's construction of this normal form and shows how it answers 
the two questions above, for a pair of coupled oscillators with linear part as in (1.5). 
A general exposition of these ideas of Arnold is given in [43, pages 305-320]. 

The reduction to Jordan form of any matrix A is accomplished by a similar- 
ity transformation. Consider the Lie group GL(n,C) of all nonsingular matrices 
SeC . Then the group of similarity transformations on C™ has the elements 
Ad s , S £ GL(n,C), defined by 

(2.2) Ad s {A) = SAS~\ A 6 C™ 2 , S 6 GL(n, C) . 
The group orbit of a given A is 

(2.3) 0{A) = {Ad s {A) = SAS- 1 | 5eGL(n,Q} 

2 

and this is a smooth submanifold of C™ . The orbit of A thus consists of all matrices 
similar to A. It is well known that similar matrices share the same eigenvalues, 
with the same multiplicities and the same Jordan blocks (perhaps permuted). A 
convenient orbit representative for any given A is its Jordan form, with Jordan 
blocks ordered in decreasing size, for example A c in (2.1). 

Now, combining the concepts of deformations and similarity, we say that two 
deformations B(ji) and A(ji) of the same matrix 5(0) = A(0) = Ao are equivalent 
if there exists a deformation C(//) of the identity matrix, for /j, £ M, such that 

(2.4) B( l ,) = C(,i)A( l ,)C( l i)- 1 . 

Next, this equivalence is extended to allow two different parametrizations for B 
and A. Let M and M be neighbourhoods of the origin in parameter spaces C k and 
C', respectively. Let ip : M — > AT be holomorphic at 0, and <p(0) = 0. Let A(a) be 
a family of matrices on Af . Then the family (deformation) induced from A by ip is 
the family (deformation) <p o A 

{poA){n)=A{p{n)). 

Definition. A deformation A(a) of Ao is called versal if every deformation B(ji) 
of Ao is equivalent to a deformation induced from A(a): 

(2.5) B(^ = C(^A(p(^)C- 1 (^, C(0) = /, p(0) = 0. 

A versal deformation is universal if mapping p> is determined uniquely by deforma- 
tion B(ji). It is minwersal if the dimension of the parameter space Af is minimal 
among all versal deformations. This minimum dimension is called the codimen- 
sion of the matrix Aq. The parameters in a miniversal deformation are also called 
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unfolding parameters, and a (uni-, mini-)versal deformation is sometimes called a 
( uni-, mint- jversal unfolding. 

It is clear that such a versal deformation exists and an upper bound on the 
codimension is n 2 , since one could take as parameters all the elements of an n x n 
matrix. A lower bound is given by n, for the case of distinct eigenvalues and 
a diagonal Jordan form. In general, we need a manifold transverse to the group 
orbits. To this end, compute tangents to the group orbits in (2.3), for fixed A = Ao- 
The derivative of Ads(ylo) = SA^S -1 , with respect to S at the identity I, for fixed 
Aq, is easily found to be the linear operator defined by commutation with Aq 



(2.6) 



DAdj : C 1 C , DAdjS = [S, A ] = SA - A S. 



The centralizer Za of matrix A 6 C™ is the set of all matrices which commute 
with A 

Z A = {B | [B, A] = BA — AB = 0} . 

Thus, the nullspace of DAd/ is the centralizer of Aq. The rank of DAd/ is the 
dimension of the group orbit of Aq. Since DAd/ is square, the dimension of the 
centralizer is the codimension of the group orbit, which is the minimum number of 
parameters in a versal deformation. This codimension was given by Arnold [6] in 
the following formula 



(2.7) 



EE(2i-l)n.j(A i 



= 1 3 = 1 



l(A,- 



are distinct eigenvalues of Ao , U is the number of Jordan 
> • • • > nn t (ai) are the corresponding orders of the Jordan 



where A 8 - , i = 1 , . 
blocks of Xi and n 
blocks for each A 8 -. 

Furthermore, with respect to the Hermitian scalar product (A, B) = tr(AB*), 
the orthogonal complement of the tangent to the orbit of Ao in the tangent space 
C™ is given by {B* \ B £ Z^ }. Here B* denotes the conjugate transpose of B. 
This is the form that the Fredholm Alternative takes in this setting. From this, 
Arnold obtains explicitly the versal deformations, see [6, 7, 9, 43]. 

Applying Arnold's method to the Jordan matrix A c above gives the centralizer 
Za c = {C(a)} and versal deformation A c (a 

/ «2 «1 
«2 






(2.8) 



C(a) 



explicitly as the families: 
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(2.9) 



Here a 



I «2 






1 









i + a 3 




J 



1, 



—i + «4 J 

, 4 are arbitrary complex parameters. This versal deformation 



is not unique. Rather than the above orthogonal complement to the tangent space, 
we could choose any transversal subspace of minimum dimension, and obtain a 
miniversal deformation. One convenient choice replaces the «2 entry in the upper 
left corner of A c (a) by 0, and thus reduces the number of nonzero entries by one. 
We adopt this choice. 
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Furthermore, since the original system is real, it is desirable to decomplexify 
(a). The result of decomplexification (for this choice) is the real family 



(2.10) 



/ 





1 





\ 




Cti 


«2 
















«3 


— 1 — «4 


V 








1 + «4 


«3 / 



A (a) 



, 4 have been redefined as real parameters. We call the family 



where ay, j = 1, 

in equation (2.10) the real Arnold normal form of the original matrix Ao. It has 
the desired versality property, that any deformation B(ji) of Ao is equivalent to a 
deformation induced from the family A a . In other words, any such B(ji) can be 
obtained by means of holomorphic mappings of parameters combined with similar- 
ity transformations, as given in (2.5). Among all such versal deformations, it has 
the minimum number of parameters (i.e. is miniversal). 

Finally, recall that our matrix is the linear part of a differential equation, in 
which a rescaling of time has the effect of multiplying the matrix by a positive 
real number. We choose to rescale time to keep the imaginary part of the complex 
eigenvalue equal to ±i. This has the effect of making «4 = in (2.10) and rescaling 
the other parameters. Therefore, in this context the codimension is 3, and the 
versal deformation is 



/ 1 



(2.11) 



A (a) 








«2 









1 






-1 

«3 



where now ay, j 



1,2,3 are redefined real unfolding parameters. 
It remains to establish a connection between the unfolding parameters ay in 
Arnold's versal deformation (or normal form) (2.11) and the physical parameters in 
the original system (1.5). Since the physical parameters can enter only in certain 
ways, determined by the mechanical structure of the system, it is not clear that 
the notion of versality for general differential equations is appropriate for this class. 
If in fact the physical parameters can not be transformed to the Arnold versal 
deformation parameters, then two possibilities must be considered: either there is 
a constraint on the physical system which should be made explicit, or an allowable 
parameter has been overlooked. 

For the present example, it is easy to find an explicit relationship between the 
physical and unfolding parameters. This may be done by comparing coefficients of 
the characteristic polynomials. The relationship is found to be 

-£}+■■■ 



(2.12) 



cti 

«2 
«3 



Si 



\s 2 + 



where • • • represents higher order terms. Since this mapping is invertible at the 
origin, we have proven the following. 

Proposition. The linear part of the physical model in (1.5), with parameters 
{£1,81,82}, is a miniversal deformation in the sense of Arnold. The three physical 
parameters are related to the unfolding parameters in (2.11) by equations (2.12). 

This justifies the choice of three physical parameters (£1,81,82) in equation 
(1.3), in the following sense. Although there may be other physical parameters 
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which affect the behaviour of the model, those other parameters can not lead to 
any new local behaviour which is not already included in the universal deformation 
(2.11). Thus, the 3-parameter family of matrices (2.11) is the appropriate starting 
point for the study of phenomena associated with coupled oscillators near : 1 
resonance. As observed by Arnold in [9, page 221]: 

. . . the natural object of study is not the individual object (say, 
a vector field with a complicated singular point), but a family so 
large that the singularity of the type under consideration does not 
disappear under a small deformation of the family. This simple 
argument of Poincare shows the futility of such a large number of 
studies in the theory of differential equations and in other areas of 
analysis that it is always somewhat dangerous to mention it. 

3. Poincare Normal Form 

The Poincare normal form may be thought of as an extension of the ideas of the 
previous section from the linear case to the case of nonlinear vector fields. Arnold 
gives an excellent account of the theory and computation of the Poincare normal 
form in [9, Chapter 5]. Other useful references are [10, 12, 14, 18, 20, 22, 23, 
28, 43]. 

The basic idea of Poincare is to perform a sequence of near-identity transfor- 
mations of the form 

(3.1) x = y + h fe (y), A = 2, 3,... 

where is a homogeneous polynomial of degree k. For each k, is chosen to 
eliminate as many as possible of the terms of degree k in the vector field. Terms 
which remain lie in a complement to the tangent space, at the identity, of the 
manifold of group orbits corresponding to (3.1), just as in the linear case of Sec- 
tion 2. The Poincare normal form is not unique, since it depends on this choice of 
complement. Furthermore, the resulting power series does not converge in general; 
nor does the sequence of near-identity transformations (3.1) converge. This lack of 
convergence is often not a problem, since in many cases the phenomena of interest 
are determined by a finite truncation of the Poincare normal form, of low degree. 

For the system of differential equations in this paper as given by (1.7) with 
H = , the Poincare normal form was first computed by Moson [31, 32]; see also 
[10, 12]. The Poincare normal form inherits the symmetry property (1.2), and can 
be chosen to have an additional normal form symmetry which is the S 1 symmetry 
inherent in the Andronov-Hopf bifurcation [20]. The resulting normal form may 
be written (formally) as 

f yi\ ( l \ / j/i \ / fa yi 

(3.2) [ y 2 = 2/ 2 + 4>m + foy2 
\ z ) V i J \ z J + i<f> 4 ] z 

to which may be added the conjugate of the last equation, for z. Here <f>j, j = 
1, . . . ,4, are real valued formal power series (which may not converge) in the Hilbert 
basis of invariant polynomials {y\, zz = y\-\-y\ }, i.e. 

<t>3 = ^jiy'l,^), i = i, • • • ,4. 

Notice the formal resemblance between (3.2) and the complex Arnold normal form 
in the linear case, given by (2.9). In fact, if we were to allow the terms involving <f>j 



10 



W.F. LANGFORD AND K. ZHAN 



in the vector on the right-hand side of (3.2) to have linear parts, then the linear part 
of the Poincare normal form would be the Arnold normal form of its linearization. 

The Poincare normal form (3.2) can be further simplified by the near-identity 
transformation: 

(2/1, 2/2, z) (2/1,2/2 + hVi,z), 
which eliminates the nonlinear terms in the 1)1 equation of (3.2). Now, the Poincare 
normal form matches the Arnold normal form as given in (2.10), in the sense that 
we can identify 

^• = ^(0,0), i = l,...,4. 

Combining in one equation the versal deformation of the linearization, i.e. the 
Arnold normal form given by (2.11), and this simplified Poincare normal form, 
yields 

/ in \ / 1 \ / yi \ I 
(3.3) 2)2 = «i «2 2/2 + ^12/1+^22/2 

\ z ) \ i + a 3 J \ z J \ [<j> 3 + i<j> A ]z 

where 2/1, 2/2 £ 1, z £ C, and <f>j = <f>j(y 2 , zz), <f)j(0,0) = 0, j = 1,... ,4, and we 
have rescaled time to make ia^ = 0. 

Finally, we decomplexify the z-equation to real form, but in polar rather than 
cartesian coordinates, to better exploit the normal form symmetry. Let 

z = re is . 

Then (3.3) can be written 

2/2 

[«i + <f>i]yi + [a 2 + <f>2]y2 
[«3 + $3] r 
l + <f>4 ■ 

Note that 6 does not appear in the right hand sides of equations (3.4), so the 6- 
equation is decoupled from the other three. This is due to the S 1 symmetry of the 
normal form, induced by the Andronov-Hopf bifurcation in z, and it effectively re- 
duces the problem from four dimensions to three. Given a solution (t/i (t), 2/2^), r (^)) 
of the first three equations of (3.4), then the solution 6(t) is given by 

6{t) = 6 +t+ f Myi(s),r 2 (s))ds . 
Jo 

All but exotic behaviour of this system is already determined by a cubic truncation 
of the three-dimensional Poincare normal form. This can be written 

2/1 = 2/2 

(3.5) 2)2 = [«i + <j>i, 102/1 + ^i,oiJ* 2 ]2/i + [«2 + (t>2, 102/1 + <t>2, oi?* 2 ]2/2 

r = [a 3 + (fe, 102/1 + <f>3,oir 2 ]r 

where 4>j t ki is the coefficient of the term y'f k r 21 , in the power series expansion of 
4>j(y 2 , r 2 ). The 6 equation has been omitted in (3.5), but has a similar expansion. 
The system (3.5) is explored in the next section. 

By standard arguments, hyperbolic solutions of (3.5)are preserved as solutions 
of (3.4), and in fact correspond to higher dimensional solutions of the original four 
dimensional system. For example, hyperbolic equilibria of (3.5) correspond to hy- 
perbolic periodic orbits of the original system. However, as has been shown in 



(3.4) 



2/1 

2/2 

r 
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[18, 25, 26, 43], in special cases the higher order remainder or "tail" of the nor- 
mal form can dramatically change the dynamics from the behaviour of a truncated 
normal form such as (3.5), since the tail breaks the normal form symmetry which 
led to the decoupling of the 6 equation. For example, when (3.5) has homoclinic 
or heteroclinic orbits, then the tail effects can produce homoclinic tangles, Smale 
horseshoes and chaotic dynamics. The singular nature of these normal form pertur- 
bations means that splitting of separatrices will be associated with exponentially 
small splitting problems, see [38, 19, 40]. These considerations will be left for 
investigation elsewhere. 

4. Bifurcations of Codimension One and Two 

Identifying the collection of all possible phase portraits of this system near 
the codimension three bifurcation point, in the sense of all equivalence classes with 
respect to topological equivalence, remains an open problem. However, the stratified 
subvariety (as defined in [9, page 227]) of bifurcations from the trivial solution in 
the parameter space M 3 takes a remarkably simple form, and is presented in this 
section. 

The stability of the trivial equilibrium solution is determined by the eigenvalues 
of the linear part Ao(a), which are found to be 

OL2 I ( OL2 \ 2 

(4.1) Ai >2 = — ± \jai + J , A 3)4 = a 3 ±i. 

All four eigenvalues have negative real part if and only if ay < 0, j = 1,2,3. 
Therefore, the trivial equilibrium solution is asymptotically stable, strictly inside 
the negative octant of the parameter space M C M 3 , see Figure 1. As the point 
(«i, «2, 0:3) leaves the negative octant, bifurcations occur. These bifurcations from 
the trivial equilibrium solution are called primary bifurcations. There are seven 
possibilities, see Figure 1: a point may leave by crossing through one of the three 
coordinate planes (codimension one), through one of the three coordinate axes 
(codimension two) or through the origin (codimension three). A generic exit is 
through a plane, codimension one. However, interesting and important mode- 
interactions occur around the intersections of the planes. These are best understood 
by studying the codimension two bifurcations on the coordinate axes. All six of 
these codimension one and two bifurcations have been studied previously, and are 
well understood in general. Their manifestations in the present system are described 
in the following subsections. 

4.1. Primary Bifurcations of Codimension One. The three primary bi- 
furcations of codimension one are: an Andronov-Hopf bifurcation with frequency of 
order 0(1) which we label Hi, an Andronov-Hopf bifurcation with small frequency 
0(|ai|), labeled Ho, and a pitchfork bifurcation labeled P; see Figure 1. Recall 
that the unfolding parameters in Figure 1 are related to the physical parameters 
H = (£1,81,82) by equations (2.12). 

4.1.1. Primary Andronov-Hopf Bifurcation Hi. For the system (3.3), the z— 
plane is a two-dimensional invariant subspace, on which the normal form equations 
reduce to 

2/1 = 2/2 = 

(4.2) r = [a 3 + (f> 3 (0,r 2 )] r 

6 = l + ^ 4 (0,r 2 ). 
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Figure 1. Stratified subvariety of primary bifurcations. 



This has the family of periodic solutions 

2/1 = 2/2 = 

( 4 -3) r = + - 

9{t) = 6» + [1 + ^(0, r 2 )]t (mod 2tt) . 

This is a classical Andronov-Hopf bifurcation, giving birth to a limit cycle in the 
(t/3, t/4)-plane with period near 2ir (angular frequency 1) as the parameter 03 crosses 
through 0, provided 0:1,0:2 7^ 0. Generically ^3^1 7^ and the bifurcation is 
nondegenerate, which is assumed to be the case in this paper; then the limit cycle 
is locally unique. Degenerate cases of Andronov-Hopf bifurcation are explored in 
detail in [17], but will not be considered here. If ^3^1 < 0, then the bifurcation 
is supercritical (the limit cycle exists for 03 > 0) and if ^3^1 > the bifurcation 
is subcritical (exists for 03 < 0). A supercritical limit cycle is stable if 01,02 < 0, 
otherwise the cycle is unstable. 

4.1.2. Pitchfork Bifurcation P. Any equilibrium solution of (3.3) lies on the y\ 
axis, and satisfies 



(4.4) 



[ai + <£i(s/?,0)]s/i =0 
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This is a classic pitchfork bifurcation, yielding in addition to the trivial solution 
y = 0, a symmetric pair of nontrivial solutions 



branching from as ct\ crosses through zero, provided ^1,10 7= and «2,«3 7= 0, 
which we assume. The pitchfork is supercritical if 0i,io < and subcritical if 
^1,10 > 0. Supercritical solutions are stable if «2,«3 < 0, otherwise both are 
unstable. 

4.1.3. Primary Andronov-Hopf Bifurcation Ho. The Poincare normalform (3.3) 
leaves the (t/i, t/2)~plane invariant. On this plane, the normal form reduces to 

Vi = 2/2 

(4.5) y 2 = [ai + <Ms/?,0)]s/i + [a 2 + <M!/i,0)]s/2 

z = . 

For fixed ct\ < 0, this system has an Andronov-Hopf bifurcation as «2 crosses 
through zero, for any a 3 ^ 0. The new limit cycle has a small frequency, asymptotic 
to y/— ct\ . The change of variables x = yi, y = —y'j/oj, where oj 2 = —ct\ transforms 
(4.5) to the planar system 



i\_(0-u\(x\ ( 



This planar system has the polar normal form, with x + iy = pe llp 
(4.7) p = 



" 2 , 1 J, 2 
Y + g02,lOP 



<p = y-«i - \y ) +b p +■■■ ■ 

From this it follows that the limit cycle has amplitude given by 



/ 4« 2 | 

Y ^2,10 

and is super-(sub-)critical if 02,io < (> 0). 

4.2. Primary Bifurcations of Codimension Two. The Poincare normal 
form (3.3) has three bifurcations of codimension two, from the trivial solution in 
state space. In the parameter space these occur along the three coordinate axes, 
see Figure 1. 

4.2.1. Bogdanov-Takens Bifurcation BT. As observed for equation (4.5), the 
(yi ; £/2)~plane is invariant for the flow of the normal form. On this plane, for 
a 3 zfz and («i, a^) in a neighbourhood of (0, 0), the system has a Z^-symmetric 
Bogdanov-Takens bifurcation [11, 41], see line BT in Figure 1. It yields stable solu- 
tions only if «3 < 0. In addition to the codimension one primary bifurcations P and 
Ho described above, there are rich possibilities including simultaneous secondary 
Andronov-Hopf bifurcations from the pitchfork pair of equilibria, coalescence of 
two limit cycles, orbits homoclinic to the origin and heteroclinic orbits joining the 
two pitchfork equilibria. A detailed description of the versal deformations, bifurca- 
tion diagrams and phase portraits is given by Arnold [9, page 300], and need not 
be repeated here. 
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New possibilities arise in this system via the interaction of the Bogdanov- 
Takens bifurcation with the primary Andronov-Hopf bifurcation Hi , leading to a 
flow in a phase space which reduces to M 2 x S 1 . The simple homoclinic orbits of the 
planar Bogdanov-Takens bifurcation then lead to transversally intersecting two- 
dimensional stable and unstable manifolds of periodic orbits and thus to chaotic 
dynamics, as for example in [18, 28, 42, 44]. These phenomena are under further 
investigation. 

4.2.2. Pitchfork-Hopf Mode Interaction PH. On the 02-axis (0:1,0:3) = (0,0), 
the system has a simple zero eigenvalue, with a Z^-symmetry on the corresponding 
eigenspace t/i, and a pair of purely imaginary eigenvalues ±i. In a neighbourhood 
of the 02-axis for 02 < 0, the normal form (3.3) has a pitchfork-Hopf mode inter- 
action, see PH in Figure 1. One finds the expected pitchfork and Andronov-Hopf 
primary bifurcations P and Hi, and secondary bifurcations of mixed-mode solu- 
tions involving nonlinear coupling of these primary bifurcations. A complete listing 
of these primary and secondary bifurcation branches for nondegenerate pitchfork- 
Hopf mode interactions is given in [27], see also [18]. There is a possibility also 
of tertiary bifurcation to a quasiperiodic flow on an invariant torus with a small 
second frequency [25, 39], discussed further below. 

To investigate these mode-interactions, in (3.4) we set y\ = 1/2 = 0, r = 0, 
suppress the 6 equation and retain only leading order terms in <f>i and ^3. Then 
equation (3.5) reduces to 



Note that the nonzero quantities on the right must be positive for the solutions 
to be real. The first three of these solutions are the previously identified trivial, 
primary pitchfork P and primary Hopf Hi solutions, respectively. The fourth is 
a mixed-mode solution, which joins the P and Hi solutions, at the two ends of 
its interval of existence, in secondary bifurcations. Due to symmetry, the primary 
bifurcations have two symmetric branches ±2/1, ±r, and the secondary bifurcations 
have four branches; however, in the full dynamics the solutions with ±r give the 
same periodic orbit with a phase shift. The secondary bifurcation from P is an 
Andronov-Hopf bifurcation, while the secondary bifurcation from Hi is a pitchfork 
bifurcation of two new limit cycles. While all of these results were deduced only from 
the cubic truncation (4.8), they extend by standard arguments to the full equations, 
using assumptions (4.9), the Implicit Function Theorem and transversality. 



(4.8) 



= [oi + 1,1 02/1 + <^i,oij* 2 ] yi 

= [a 3 + <fe, 102/1 + <f>3,oir] r . 



For nondegeneracy, we assume 

(4.9) Ai = ^1,10^3,01 - ^3,10^1,01 ^ 0, ^1,10^3,01 ^ . 

Then (4.8) has four branches of solutions 



(j/ 2 ,r 2 ) = (0,0) 
(4.10) (y\y) = 
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When higher order terms are included, in the case ^1,10^3,01 < it is possible 
to find a tertiary bifurcation, from the mixed-mode secondary branch above to a 
flow on a two-dimensional invariant torus with a small second frequency [13, 18, 
25, 39, 42]. Explicit conditions for the existence of this type of 2-torus bifurcating 
from a steady state - Hopf mode interaction, and a description of the flow on the 
torus, have been given by Scheurle and Marsden, see [39]. It should be noted that 
this 2-torus exists only very locally in parameter space, in contrast to the 2-torus of 
the next subsection. The secondary radius of this 2-torus grows very rapidly from 
bifurcation, and it becomes sensitive to the effects of the symmetry-breaking tail 
(or remainder) of the normal form, which may lead to chaos. 

4.2.3. Hopf-Hopf Mode Interaction HH. In a neighbourhood of the negative ct\ 
semi-axis (see HH in Figure 1), nonresonant Hopf-Hopf mode interactions occur, 
involving the Hi and Ho bifurcations. For (0:2, 0:3) = (0, 0), the full system has two 
pairs of purely imaginary eigenvalues, which are nonresonant because the frequency 
ratio is very small. An introduction to the dynamics of nonresonant Hopf-Hopf 
mode interactions is given in [18, Chapter 7]. The four-dimensional system can be 
written as two pairs of amplitude-phase equations, in (r, 9) and (p, <p), corresponding 
to the two primary Andronov-Hopf bifurcations Hi and Ho respectively. The two 
phase equations decouple due to nonresonance, leaving two amplitude equations 
which have the TL'i © TL'i symmetric normal form 



(4.11) 



r 
P 



a 3 

i«2 



+ cp 2 + <f> 3 ,oir 2 + ■ ■ ■] 



\h,wp' 2 + dr' 2 + ■ ■ ■} 



This has four families of solutions satisfying (r,p) = (0,0), given, subject to the 
obvious nondegeneracy conditions, to leading order by 



(4.12) 



(r, P) 
{r, P) 

{r, P) 
{r, P) 



(0,0) 




where A2 = g</ ) 2,io</ ) 3,oi — erf. The first three of these are the previously obtained 
trivial, Andronov-Hopf Hi and Andronov-Hopf Ho solutions, respectively. The 
fourth gives a new family of two-frequency solutions, on an invariant 2-torus. The 
local persistence of this invariant torus is guaranteed by standard arguments in- 
volving normal hyperbolicity and nonresonance. This 2-torus differs from the 2-tori 
discussed in the previous subsection, in that it is more robust, and exists over a 
much wider region of parameter space. Thus, one may expect that this torus may 
be observed more easily in physical systems. 

This 2-torus appears through a secondary bifurcation from either of the two 
primary Andronov-Hopf families, and may smoothly link these two primary fami- 
lies. Generically, this family of 2-tori is either always stable or always unstable. To 
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leading order, the 2-torus bifurcates from the Hi periodic orbit when 



1 I Q/3 
(4.13) a 3 d = -a 2 <f> 3 ,01, P = 0, r = . /-- . 

2 V 03, 01 

and bifurcates from the Hq periodic orbit when 



1 / 4«2 

(4.14) a 2 c = -a 3 02,io, P=\--; , r = 0. 

4 y 02,10 

For parameters near (4.13), the solutions on the 2-torus are close to the Hi periodic 
orbit with primary frequency 0(1), but have a small amplitude modulation with 
frequency 0{\/—ai). The amplitude of the modulation grows along the branch 
from (4.13) towards (4.14). This behaviour resembles that observed by Nayfeh et 
al. [34, 35, 36]. In certain cases, there is a further possibility of a stable 3-torus 
bifurcating locally from the 2-torus [21], and of chaotic dynamics; however, the 
analysis of these phenomena requires consideration of higher order terms and will 
not be pursued here. 

5. Example 

The analysis of the previous section makes it easy to find an example of a 
pair of coupled oscillators, of the general form (1.3), which exhibits (for a range of 
parameter values) the behaviour described by Nayfeh et al. [34, 35, 36]. One such 
example is 

(5.1) x + Six + i2x + (x + x)(x 2 + y 2 ) = 

y + S 2 y + y+{y + y){-0.2x 2 + y 2 ) = 

with parameter values chosen as SI = 0.3, Si = —0.2 and S 2 = —0.25. Equations 

(5.1) may be interpreted as a pair of nonlinearly coupled Duffing - van der Pol 
oscillators, as follows. In the first equation, the x 3 term is the Duffing nonlinear 
restoring force and xx 2 is the van der Pol nonlinear damping force. These are both 
modified through nonlinear coupling effects of the second oscillator, represented 
by the xy 2 and xy 2 terms, respectively. An analogous interpretation holds for the 
second equation. 

In accord with the analysis of the previous section, solutions of this system 
approach an invariant 2-torus, which exhibits a "fast oscillation" in the (y, y) vari- 
ables together with a "slow modulation" in the x variable; see Figure 2 for two 
such solutions, which have been computed numerically and plotted using Maple. 
The physical parameters chosen in this example correspond to values of the un- 
folding parameters near the negative «i-axis in Figure 1, given by equation (2.12) 
approximately as 

(5.2) ai = -0.3, a 2 = 0.2, a 3 = 0.125 . 

A more detailed study of this phenomenon, for more physically relevant examples, 
and including both numerical and analytical results, is in progress. 

6. Conclusions 

We may draw several conclusions from this study. The first important conclu- 
sion is that the simultaneous occurrence of Andronov-Hopf and Bogdanov-Takens 
bifurcations is a codimension-three phenomenon, requiring three independent pa- 
rameters for its exploration and understanding. In the versal deformation of the 
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Figure 2. Two solutions approaching the invariant torus of Sec- 
tion 5, from two initial points inside and outside the torus, 
(x(0),i(0),y(0),y(0)) = (0.1,0,0.1,0) and (-0.3,0,0.4,0.9), re- 
spectively. 

codimension three bifurcation, there are primary bifurcations of six different types; 
three of codimension one and three of codimension two, all of which are classi- 
cal. The stratified subvariety of primary bifurcations is shown in Figure 1. The 
three parameters in Arnold's versal deformation may be chosen so that they di- 
rectly determine each of the three primary bifurcations of codimension one: the 
plane a\ = determines the primary pitchfork bifurcation P, a 2 = determines 
the primary Andronov-Hopf bifurcation Ho with low frequency (of order \/—a\) 
and «3 = determines the primary Andronov-Hopf bifurcation Hi with high fre- 
quency (of order 1). The three codimension two bifurcations correspond to the 
three coordinate axes in Figure 1. They lead to secondary bifurcations of new 
families of solutions which are linkages between the primary bifurcation families of 
codimension one, and to further rich and complex dynamics. 

The three versal unfolding parameters are related by a nonsingular transforma- 
tion to the three physical parameters of the coupled oscillators in the Introduction, 
by formulae (2.12). This verifies that a versal deformation gives the proper setting 
for investigation of the possible behaviours of the class of coupled oscillators defined 
in the Introduction. 

The nonlinear coupling between the primary high frequency and low frequency 
Andronov-Hopf bifurcations can lead to a quasiperiodic flow on a 2-torus, with its 
two frequencies close to those of the two primary Andronov-Hopf bifurcations. This 
gives a plausible explanation for the observations of Nayfeh et al. [34, 35, 36], of 
low frequency modulation of a high frequency mode, and of energy transfer from 
high to low frequency modes, near : 1 resonance. The case studied by Nayfeh 
corresponds to a choice of unfolding parameters a\ < (fixed) and «2,«3 both 
small and negative; i.e., near the negative a\ axis marked HH in Figure 1. This 
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is the region where existence of a 2-torus is predicted by the normal form analysis 
and confirmed by the numerical analysis in Section 5. 

The present study yields several indicators of directions for future research. 
Systems of coupled oscillators with less symmetry than required by (1.2) are im- 
portant in applications, and may be investigated by the same normal form methods 
as used here. The effect of higher order terms, particularly in the onset of more com- 
plex dynamics, is worthy of investigation. For example, in the standard symmetric 
Bogdanov-Takens analysis, the double zero eigenvalue may lead to homoclinic and 
heteroclinic orbits in the phase plane. When coupled with the primary Andronov- 
Hopf bifurcation Hi as here, it will lead to phenomena of homoclinic chaos, which 
merit further investigation. A more detailed investigation of the mechanism of 
energy transfer from high to low frequency modes, for more realistic physical ex- 
amples, and using both analytical and numerical techniques, is in progress. As 
well, the application of these results to concrete systems in engineering, biology 
and other fields may lead to new insights. 
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